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ABSTRACT 


A constellation  of  satellites  that  gives  multiple  con- 
tinuous coverage  of  the  whole  earth  is  studied  here.  A com- 
puter program  is  described  which  helps  to  find  the  optimal 
configuration.  In  particular,  performance  of  a twenty- four 
satellite  constellation  is  analyzed  for  triple  coverage  and 
the  results  given.  Constellations  where  the  number  of  orbital 
planes  equal  the  number  of  satellites  (N  = P)  and  constellations 
where  N ^ P are  both  analyzed.  Circular  orbits  are  assumed. 

The  computer  program  is  relatively  short  using  a sampling  grid 
over  a quarter  of  the  globe.  The  results  are  approximately 
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Introduction 


In  order  to  update  a vehicle's  state  vector,  single, 
double,  triple  or  quadruple  coverage  may  be  desired  from  a 
constellation  of  artificial  satellites.  One  important  con- 
sideration to  the  constellation  designers  is  the  optimal 
configuration  in  terms  of  pure  geometrical  distribution. 
Thus,  what  is  the  best  coverage  N satellites  can  produce 
and  what  are  the  Keplerian  elements  of  each  satellite  at 
some  point  in  time? 

In  the  literature  one  can  find  the  optimal  coverage 
for  constellations  of  up  to  15  satellites^-.  There  is  pres- 
ently a need  for  better  coverage  than  is  afforded  by  these. 
The  earth  central  angle  between  user  and  satellite  is  too 
large  for  many  applications. 

J.  G.  Walker  has  laid  down  much  of  the  ground  work  for 
analysis  of  such  constellations.  What  is  needed,  however, 
is  a method,  employing  a computer  program,  to  simplify  and 
make  possible  analysis  of  systems  where  N is  considerably 
greater  than  15. 

Constellations  where  each  satellite  is  placed  in  a 
single  orbital  plane  (N  = P)  and  those  where  equal  numbers 
of  satellites  are  contained  in  fewer  planes  (N  > P)  are  both 
considered  here.  For  the  N = P case,  type  I,  the  inclina- 
tion is  assumed  the  same  for  all  orbits,  the  spacing  in 


longitude  of  ascending  node  (LAN)  is  made  equal  to  360°/N 
and  the  phase  angle  between  satellites  is  given  by  F x 360/N. 
F is  an  integer  whose  optimal  value,  along  with  inclination, 
is  to  be  found. 

For  the  N > P case,  type  II,  all  possibilities  where  P 
is  a factor  of  N are  considered.  The  increment  of  LAN  is 
then  360/P,  the  phase  angle  is  defined  the  same  way  as  for 
type  I constellations,  and  the  inclinations  for  all  P planes 
is  the  same. 

The  'Grader'  and  'Grader  1'  programs  analyze  these  two 
types  of  constellations  respectively.  They  take  a grid  of 
samples  over  half  of  the  northern  hemisphere  to  give  an  ap- 
proximation of  the  performance  of  the  different  possible 
constellations.  It  is  meant  to  be  an  aid  to  the  systems 
engineer  who  wants  to  get  an  idea  of  what  the  constellation 
will  look  like.  One  thus  arrives  at  the  most  probable  opti- 
mal configuration,  the  number  of  orbital  plane  (P) , the 
phase  angle  number  (F) , and  the  largest  earth  central  angle 

(R  ) between  the  user  and  the  satellite  which  must  be  tol- 
max 

erated  to  get  the  coverage  specified. 


Section  I 


A Look  at  Whole  Earth 
Continuous  Coverage 

In  'whole  earth  continuous  coverage'  one  is  striving  to 
have  at  any  point  in  time,  anywhere  on  the  globe,  a prescribed 
number  of  satellites  'visible'.  The  'prescribed  number'  de- 
pends of  course  on  the  user's  intent.  If  one  desires  single 
coverage  then  at  least  one  satellite  should  be  visible,  and 
so  on.  Visibility  requirements  can  usually  be  very  neatly 
defined  by  a maximum  allowable  earth  central  angle  between 
the  user  and  the  satellite.  In  the  example  chosen  in  this 
paper,  we  are  looking  for  triple  coverage  with  a maximum 
earth  central  angle  of  approximately  45  degrees. 

Note  that  although  the  coverage  requirement  may  be  ful- 
filled, the  type  of  coverage  the  satellites  afford  could  be 
poor  for  navigation.  A measure  of  the  navigability  of  a 
given  constellation  is  not  included  here  but  is  being  studied 
at  this  time  ^ . 

Results  of  the  study  are  shown  in  Table  1.  The  N/P/F 

notation  denotes  the  number  of  satellites  (N)  , tile  number  of 

orbital  planes  (P) , and  the  phase  angle  number  (F) . R 

rndx 

represents  the  largest  angular  distance  to  acquire  three 

satellites.  This  number  should  be  as  low  as  possible.  The 

inclination  (6)  is  to  a reference  plane  which  can  be,  but 

does  not  have  to  be,  the  equatorial  Diane.  D . represents 

min 
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Table  1 
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Summary  of  Triple  Coverage  Analysis 


Constellation 

Code 

Minimum 

R 

max 

At  Incli- 
nation 6 

Minimum  Distance 

Between  Any  Two 

Satellites  (D_.  ) 
mm 

(deg) 

N/P/F 

(deg) 

(deg) 

24/24/14 

51 

53.0 

17.1 

24/3/1 

50 

53.3 

12.2 

the  minimum  angular  distance  between  two  satellites  encoun- 
tered with  this  scheme.  Constellations  where  D . was  equal 

mm 

to  zero  were  discarded  due  to  the  relatively  high  probability 
of  collision  for  circular  orbits  of  equal  altitudes.  Poor 
navigation  also  results  at  those  instances  and  places  where 
this  occurs.  Ualker  also  takes  this  approach.  The  assump- 
tion is  made  that  the  proper  station  keeping  is  available  to 
keep  the  relative  geometry  constant.  Similar  perturbations 
to  all  the  orbits  may  be  ignored  such  as  those  shifting  the 
longitudes  of  ascending  nodes. 
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Section  II 

Finding  the  Optimal  Configuration 

Finding  the  optimal  configuration  can  be  long  and  tedious 

if  attempted  by  purely  deterministic  approaches,  especially 

when  the  number  of  satellites  is  large.  There  are  too  many 

variables  involved  as  can  be  supported  by  anyone  who  has  at- 
. ( 2 ) 

tempted  it  . In  order  to  speed  up  the  process,  a program 
called  'GRADER'  was  introduced  (see  Section  IV  and  Appendix  A 
for  more  details) . 

The  program  narrows  the  field  of  contenders  to  one  or 
two  configurations  for  an  N-satellite  constellation.  It 
helps  to  pinpoint  the  number  of  orbital  planes,  F number, 
and  inclination  of  the  orbits.  This  is  accomplished  by  an 
iterative  method  on  the  part  of  the  user.  It  is  not  intended 
to  give  an  extremelv  accurate  optimal  inclination  or  R 

max 

D . , the  minimum  distance  between  two  satellites,  must  also 
min 

be  found  by  other  means.  An  algorithm  which  can  be  used  in 

finding  the  actual  value  of  the  optimal  inclination  and  the 

corresponding  R„  are  outlined  in  Section  III. 
max 

The  way  the  figure  of  merit  is  given  to  a configuration 
by  the  Grader  program  will  be  discussed  in  detail  in  Section  IV 
and  Appendix  A.  Basically  the  program  looks  at  a sufficiently 
large  sampling  area,  checking  for  single,  double  or  triple 
coverage  at  each  sampling  point.  The  number  of  failures  in 
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coverage  is  divided  by  the  number  of  samples  taken  to  get 
failures  per  samples  (FPS) . To  find,  then,  the  best  perfor- 
mance possible  from  a N-satellite  constellation,  the  procedure 
shown  below  has  been  found  helpful: 

a)  One  starts  with  N = P constellations 

b)  then  runs  the  program  for  the  different  possible 
F numbers  at  elapsed  angle  equal  to  0.  The 
threshold  angle  for  visibility  must  be  set  so 
that  the  curves,  as  shown  in  Figures  1,  2,  3 

to  7,  do  not  bottom  off.  Inclinations  can  be 
set  initially  to  run  from  40  to  70  at  increments 
of  5 degrees. 

c)  The  N > P constellations  can  then  be  checked 
by  running  the  program  for  different  values 

of  P.  At  each  P,  the  meaningful  F numbers  are 
checked  also  at  elapsed  angle  set  to  zero  and 
inclination  varying  from  40  to  70. 

d)  The  configurations  with  the  lowest  minima  can 
then  be  checked  at  other  elapsed  angles  until 
the  inclination  with  the  lowest  maximum  value 
is  clearly  defined. 

For  N » 24  with  triple  coverage,  the  configuration  which 
behaves  the  best  seems  to  be  24/3/1.  This  is  found  from  fig- 
ures 1 to  7 which  show  that  although  a few  configurations 
have  very  low  FPS  values  at  E.A.  equal  to  zero,  e.g.  24/6/5, 
24/12/1,  24/24/8,  24/24/14,  24/24/16,  at  other  E.A.  values 
they  have  higher  minimums  than  *4/3/1. 
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Inclination  (deg) 

Figure  7 - FPS  vs.  Inclination 


Three  points  regarding  time  saving  measures  in  the  search 
can  be  made  at  this  point.  First,  one  does  not  need  to  vary 
the  elapsed  angle  until  all  the  F numbers  have  been  analyzed 


at  an  EA  equal  to  zero.  Only  those  F values  which  look  promis- 
ing need  to  be  investigated  further.  If  a few  values  of  F 
show  promise,  then  the  node  angle  through  which  the  satellites 
must  travel  before  the  geometry  repeats  itself  can  be  a good 
guide  for  the  next  move  in  the  search.  The  smaller  the  repe- 
tition angle,  the  better.  An  idea  of  the  repetition  angle  can 
be  had  by  looking  at  just  the  nodes  of  each  satellite  and  see- 
ing when  a satellite  will  move  to  the  equatorial  position,  etc. 
Looking  at  24/24/14,  we  see  that  with  a phase  angle  of: 

Au  * F x 360/N 
= 14  x 360/24 

= 210  deg. 


longitude 

! Of 

ascending  nodes 

and 

node  angles 

are : 

Sat 

Sat 

No. 

a 

Id 

No. 

a 

(j0 

0 

0 

0 

6 

90 

180 

1 

15 

210 

7 

105 

30 

2 

30 

60 

8 

120 

240 

3 

45 

270 

9 

135 

90 

4 

60 

120 

10 

150 

300 

5 

75 

330 

11 

165 

150 
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Sat. 

No. 

n 

u 

Sat. 

No. 

n 

oj 

12 

180 

0 

18 

270 

180 

13 

195 

210 

19 

285 

30 

14 

210 

60 

20 

300 

240 

15 

225 

270 

21 

315 

90 

16 

240 

120 

22 

330 

300 

17 

255 

330 

23 

345 

150 

Observing 

the  node  angles  of 

satellites 

5 and  17, 

at  330°, 

as  well  as 

10  and 

22  at  300°, 

and 

7 and 

19  at  30° 

, one  can 

be  fairly 

certain 

that  the  overall 

. geometry  will 

repeat 

itself  eve 

ry  30  degrees. 

Second,  as  mentioned  before. 

some 

phase  angle  numbers 

cause  satellites  to  coincide 

periodically.  The  odd  numbered 

F's  in  the 

24/24/F  constellations 

were 

found  to  have  this 

characteristic.  The  F's  versus  inclination  graph  for  24/24/5 
is  shown  in  Figure  1.  The  indication  is  that  of  high  failure 
rates  as  would  be  expected. 


Thirdly,  some  F numbers  are  unlikely  candidates  for 
optimality  judging  by  the  inherently  poor  distribution  which 
they  offer.  For  example  24/24/0  means  that  all  satellites 
will  have  the  same  node  angle,  thus  they  will  all  be  at  the 
equator  initially  and  move  upwards  at  the  same  rate.  This 
is  obviously  poor  for  global  coverage  and  need  not  be  ana- 
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lyzed  further.  F numbers  such  as  1 and  23  can  also  be  en- 
visioned as  rather  poor  but  may  still  be  checked  quickly  for 
completeness . 

Up  to  this  point  only  integer  numbers  have  been  used  for 
F.  Justification  for  this  can  be  found  in  references  2 and  3. 
They  have  found  that  fractional  values  of  F can  be  discarded. 
One  may  note  that  if  a fractional  value  were  used,  a certain 
discontinuity  would  result  in  the  distribution. 

The  different  configurations  were  initially  checked  from 
40  to  70  degrees  in  inclination  at  increments  of  5 degrees. 
This,  as  opposed  to  finer  plotting,  gives  a fairly  good  idea 
of  their  performance  while  saving  much  computer  time.  Incli- 
nations below  40  degrees  or  above  70  were  not  checked  due  to 
the  indications  of  the  curves  as  well  as  references  3 and  4. 

Figures  3 and  6 show  the  results  found  from  constella- 
tions 24/24/14  and  24/3/1  which  were  the  two  most  promising 
ones.  Note  that  for  24/24/14,  relative  to  the  sampling  grid, 
there  is  repetition  after  360  degrees  of  elapsed  angle  at 
which  time  the  satellite  numbered  0 is  back  at  the  equator. 
Therefore,  every  10  degrees  of  elapsed  angle  was  checked  up 
to  360  and  increments  of  2 degrees  were  used  in  the  area  of 
interest.  Only  elapsed  angles  of  73  to  92  are  shown  for  clar- 
ity. All  other  elapsed  angles  produced  curves  with  FPS  values 
below  those  shown. 
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Constellation  24/3/1,  which  contains  3 satellites  per 
orbital  plane,  was  only  checked  at  elapsed  angles  between 
0 and  45  degrees  since  relative  to  the  grid  pattern  this  is 
the  repetition  angle.  Note  that  this  is  confirmed  by  Fig- 
ure 6 which  shows  that  the  curves  for  E.A.  of  0 and  45  fall 
directly  on  top  of  each  other.  The  fact  that  there  are  8 
satellites  in  each  orbital  plane  at  forty  five  degrees  apart 
means  that  a satellite  at  the  equator,  for  instance,  will  be 
replaced  by  one  moving  towards  the  same  point  in  the  equator 
after  forty-five  degrees. 


Il 


The  value  of  the  largest  angular  distance  to  acquire 

three  satellites  (R  ) can  be  found  in  one  of  two  ways. 

max 

The  first  way  is  discussed  in  Section  III  and  can  be  used 
when  the  number  of  satellites  is  relatively  low  (below  24 
at  least) . It  involves  drawing  the  constellation  being  con- 
sidered on  a globe  of  some  type  and  then  finding  the  satel- 
lites surrounding  the  points  of  worst  coverage.  Using  a 

closed  form  expression  for  R , a search  is  then  made  for 

max 

the  best  inclination  and  so  forth.  The  method  is  tedious 
and  long.  It  can  be  computerized,  however  the  runs  tend  to 
be  vary  long  and  expensive  due  to  the  number  of  satellite 
combinations  which  could  be  involved  in  producing  the  worst 

4 

case.  This  has  been  done  by  Walker  . The  program,  written 

* 

in  FORTRAN,  is  approximately  700  lines  long. 

The  second  method,  employed  in  this  work,  is  to  take 


i 

\ 


A 


the  Grader  programs  and  for  the  prime  candidates,  say  con- 
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D 


stellations  24/24/14  and  24/3/1,  to  raise  the  visibility 


threshold  angle  (alpha) . Runs  are  then  made  to  determine 


what  value  of  alpha  gives  no  failures  per  sample  (FPS  ■ 0)  . 


This  aloha  value  is  a good  approximation  of  R . Shown 
■ max 


below  are  outputs  of  runs  made  for  the  constellation  24/3/1, 


Figure  8.  The  results  can  be  found  in  Table  1 of  Section  I. 


In  this  case  Rmax  is  approximately  50  degrees. 


24/  3/  1 

ELAPSED  ANGLE  = 
THREHSHOLD  AMGLc  = 


20  .00  (DEG) 
49.  (DEG) 


5 2 .99 

53.49 
53  .99 

54.49 
54.99 


FA  IL/SAMPLE 
F AIL/SAMDLE 
FA  IL/SAM^LE 
FA  I L/SAMPLE 
FA  IL/SAMPLE 


0 .COfi 
0 .00  l 1 
0.0011" 
0.0007 
C.0003 


TOTAL  NO.  CF  SAMPLES  PER  INCL.  POINT  = 2t>37 


24/  3/  1 

ELAPSED  ANGLE  = 

THPEHSHOLO  ANGLE  = 


20.00  (DEG) 
50.  (DEG) 


52.99 

53.49 
5 3.99 

54.49 

54.99 


FAIL/ SAMPLE  0.0000 
FA  I L/SAMPLE  0.0000 


FAIL/SAM=»LE  0.0000 

fail/sample  0.0000 
fail/sample  c.ooco 


TOTAL  NO.  OF  SAMPLES  PER  INCL.  POINT  = 2637 


Figure  8 - These  outputs  show  how  a threshold  angle 
of  50°  causes  all  failures  in  coverage  to 
become  zero  at  both  20  and  36  elapsed  angles, 
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24/  3/  l 

ELAPSED 

ANGLE  * 

36.00  (OEG)  <aa«*sai 

THPEHSHOLC  ANGLE  * 

49.  (DEG) 

l NC 

52.99 

FA  IL/SAMPLE 

0.00  14 

INC 

53.40 

FA  IL/SAMPLE 

0 .00  1 1 

I NC 

53.99 

FAIL/SAMPLE 

0.0007 

INC 

54  .49 

fa  il/sample 

0.0003 

I NC 

54.99 

c A IL/SAMPLE 

0.0003 

total  no.  of  sables  pep 

INCL . POINT  = 26  87 

* 

24/  3/  l 

N* 

ELAPSED 

ANGLE  = 

36.00  (DEG) 

< =3S3S3S 

THPEHSHOLC  ANGLE  = 

SO.  (OEG) 

INC 

52.99 

FA IL/SAMPLE 

c.coco 

INC. 

53  .49 

FA  IL/SAMPLE 

0 .0000 

I NC 

53.99 

FAIL/SAMPLE 

0.0000  ' 

' 

INC 

54.49 

fa  il/sample 

C.0000 

I NC 

54.99 

fail/sample 

0.0000 

TOT 

AL  NO.  OF  SAMPLES  PEP 

INCL  . PO  INT 

= 2637 

Figure  8 (continued)  Note  how  there  are  no  failures 
in  coverage  when  the  threshold  angle  is  in- 
creased from  49  to  50  degrees. 
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Section  III 


Finding  the  Point  of  Worst  Coverage 


The  method  mentioned  at  the  end  of  Section  II  and  em- 
ployed by  Walker  to  find  the  point  of  worst  coverage  will  be 
explained  below.  He  gives  a brief  explanation  of  it  in  ref- 
erence 2.  The  following  comments  might  serve  as  a help  to 
the  reader  since  the  approach  v/ill  be  covered  in  greater  de- 
tail with  some  added  modifications.  This  method  was  not  used, 
however,  in  finding  P as  the  author  had  initially  intended. 
It  was  found  to  be  somewhat  cumbersome.  In  the  end  it  does 
yield  as  precise  a value  as  is  desired.  For  the  specific 
task  at  hand  however,  an  approximate  value  is  all  that  is 
required. 

The  point  of  worst  coverage  for  a constellation  is  of 
interest  because  it  is  used  as  the  measuring  stick  for  per- 
formance. It  represents  the  greatest  angular  distance 
from  the  local  vertical  to  a third  subsatellite  point  (in 
the  case  of  triple  coverage)  anywhere  on  the  globe  and  at  any 
elapsed  angle  of  motion. 

One  might  think  that  the  point  of  worst  coverage  is  the 
one  with  the  greatest  angular  radius  to  three  equidistant 
subsatellite  points.  Looking  more  closely,  however,  one 
finds  that  moving  the  user  position  closer  to  one  of  the 
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subsatellite  points  means  that  a circle  of  radius  'R'  will 
no  longer  contain  three  subsatellite  points.  See  Figure  9. 
Thus,  the  worst  point  is  still  to  be  found.  We  are  assuming 
here  a somewhat  useful  distribution  where  the  satellites  are 
not  clustered. 


Fig.  9 The  dotted  circle  shows  worse  coverage 
than  the  dashed  circle. 

For  a true  worst  case  position,  then,  one  must  usually 

find  a point  where  three  or  more  subsatellite  points  lie  on 

the  circumference  of  a circle  with  radius  R and  two  sub- 
max 

satellite  points  lie  within  that  radius.  If  there  are  less 
than  three  subsatellite  points  on  the  circumference,  then 
there  is  a good  chance  that  this  is  not  the  worst  point. 

See  Figure  10.  Although  in  Figure  10  the  user  can  see  five 
satellites,  if  his  angular  visibility  was  only  R-e,  he  would 
only  see  two  satellites,  where  e is  much  smaller  than  R. 
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Thus  R and  not  R - e is  required  to  have  triple  coverage. 
Note  that  for  double  coverage  the  user  needs  one  subsatel- 
lite point  within  the  circle  and  so  on. 


Fig.  10  Typical  point  of  worst  coverage 
at  center  of  circle. 

For  a given  constellation,  then,  one  must  first  map 
the  satellites  (or  subsatellite  points)  on  a sphere.  See 
Appendix  B for  a program  which  gives  each  individual  satel- 
lite's latitude  and  longitude  from  the  constellation  param- 
eters N/P/F,  inclination,  and  elapsed  angular  motion.  Then 
find  the  critical  points  which  give  maximum  R values.  In 
the  case  of  24/24/14,  the  critical  points  were  though  to 
occur  at  0 and  15  degrees  elapsed  angle.  Although  this  was 
later  found  to  be  incorrect,  these  elapsed  angles  will  be 
used  as  example.  At  0 degrees,  the  satellites  producing 
maximum  R are  shown  in  Figure  11,  and  at  15  degrees,  in 
Figure  12. 
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Figure  11  - Satellites  making  up  what  was  thought 
to  be  the  worst  case  of  24/24/14  at  0° 
E.A.,  53°  inclination. 


Thus  after  the  program  'Grader'  has  placed  one  in  the 
vicinity  of  the  optimal  inclination,  a graphical  analysis  to 
find  the  satellites  making  up  the  worst  cases  is  performed. 
One  can  then  generate  a plot  such  as  in  Figure  13 . The  pro- 
gram used  for  this  is  explained  in  Appendix  C.  An  even  finer 
plot  can  be  made  to  pinpoint  the  most  favored  inclination  as 
shown  in  Figure  14  for  this  example,  that  inclination  was 
52.33°  180 


Figure  12  - Thought  to  be  worst  case  satellites. 
Elapsed  angle  is  15  degrees. 


Inclination  (deg) 


Figure  13  - A plot  of  R vs.  inclination  for  what 

were  thought  to  be  the  worst  case  elapsed 
angles  of  Constellation  24/24/14. 
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The  exception  to  the  rule  of  finding  the  worst  coverage 
is  shown  below  in  Figure  15.  In  this  case  there  are  no  sub- 
satellite ooints  within  the  circle  with  radius  R . However 

max 

there  are  enough  satellites  on  the  circumference  that  if  the 
user  moves  in  any  direction  he  still  has  triple  coverage. 


180° 


Figure  15  - Worst  case  for  24/24/8.  An  exception 
to  the  rule.  E . A.  is  60°. 

In  general,  care  must  be  taken  not  to  miss  areas  of  poor 
coverage  when  doing  this  kind  of  analysis  by  hand.  The  author 
had  access  to  a slate-covered  sphere  with  permanent  latitude 
and  longitude  lines.  This  proved  handy  indeed.  For  measur- 
ing angular  distances  an  ordinary  compass  was  used  with  chalk 
on  one  side  attached  by  means  of  rubber  bands  and  a pencil 
strapped  to  the  other  end  with  the  eraser  acting  as  the  pivot. 
Accuracies  of  +1°  or  better  were  found  when  the  graphical 
value  of  R was  compared  to  the  computed  value. 

IucLX 
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Section  IV 


The  'Grader'  Program 


The  basic  purpose  of  this  program  is  to  provide  a quick 
tool  for  checking  the  performance  of  a constellation.  It 
can  help  circumvent  hours  of  tedious  three  dimensional  plot- 
ting of  twenty  or  so  orbits  for  many  possible  phase  angles 
and  varying  inclinations.  The  program  has  been  called 
'Grader'  since  it  gives  a figure  of  merit  or  grade  to  the 
candidate  constellations. 

It  is  a simple  program  which  places  a position  vector 

in  2°  steps  in  the  latitudinal  direction.  It  stops  at  each 

latitude  and  steps  the  position  vector  in  the  longitudinal 

direction.  At  each  sample  point  it  takes  the  dot  product 

\ 

of  this  vector  and  subsatellite  position  vectors.  The  sub- 
satellite vectors  are  arrived  at  from  the  corresponding 
longitude  of  ascending  node  (LAN)  and  nodal  information  of 
each  satellite.  A transformation  is  performed  to  get  the 
information  from  this  mode  to  rectangular  coordinates  for 
the  purpose  of  taking  the  dot  product.  Dot  products,  in- 
stead of  the  corresponding  angle,  are  then  compared  to  a 
predetermined  value.  If  the  dot  product  is  greater  than 
this  threshold  value,  the  satellite  is  regarded  as  'visible' 
to  the  vehicle  (or  user) . A counter  keeps  track  of  the  num- 
ber of  visible  satellites.  Whenever  the  total  number  of 
visible  at  a particular  position  is  less  than  three,  a triple 
coverage  failure  is  noted. 
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At  the  end  of  all  the  sampling  the  total  number  of  fail- 
ures is  divided  by  the  total  number  of  samples.  Thus  one 
arrives  at  the  failure  rate  or  failures  per  samples. 

In  order  to  save  computer  time,  the  analysis  of  coverage 
is  done  over  a section  of  the  globe  and  not  over  its  entirety. 
Angles  of  90°  to  0°  latitude  and  0 to  180°  in  longitude  were 
chosen.  All  the  configurations  encountered  exhibited  symme- 
try about  the  equator  and  also  about  a plane  cutting  perpen- 
dicular to  the  equator.  This  is  the  basis  for  the  decision 
to  look  at  only  half  of  the  northern  hemisphere.  In  any  case, 
half  of  the  northern  hemisphere  is  a large  enough  area  for 
the  accuracy  desired  in  a preliminary  study  such  as  this. 

Below  is  a sample  output  of  the  Grader  for  the  24/24/14 
constellation  at  80°  elapsed  angular  motion.  It  is  plotted 
in  Figure  3. 


24/24/ 1 4 

tlAPSEO  ANGLE'  = 80. CO  (DEG)  <===== 
THPESHCLO  ANGLE  = 43.00  (DEG)  (ALPHA) 


INC 

5 C .09 

FA  I L /SAMPLE 

0.0468 

I NC 

51.99 

FAIL/SAMPLE 

0 .0446 

INC 

52  .99 

fa  IL/SAMPLE 

C .0457 

I.NC 

52L..99. 

. FAIL/SAVPLE 

C .0468 

INC 

54.99 

fail/sample 

C.0506 

The  Grader  1 program  is  almost  identical  to  Grader  only 
it  handles  cases  where  there  are  more  than  one  satellite  per 


orbit  (N  > P)  such  as  24/3/1.  The  two  programs  could  have 
been  made  into  one  rather  easily.  Shown  in  Figure  17  is  a 


sample  output  of  the  Grader  1 program. 


24/  2/  1 

EL  APS  E C ANGLE  = -•  C . -I  C < DEG  ) <=  = ====  = 

THREHSHCLD  ANGLE  = 45.  < CEG  ) 


INC 

53. ?9 

FA  IL/SAM^LS 

1.C521 

I NC 

51.9? 

c A I L/3  Ay01.  2 

•:  .143' 

INC 

52.99 

FA  !L/5A'^=>L£ 

■>.  13  37 

INC 

53.9? 

FA  XL/ 3 AM=LE 

I .'1237 

INC 

54.99 

fail/sv^u 

•:  . 3 1 

INC 

55.9? 

FA IL/SA V=LE 

:. C373 

I NC 

5o  •?  9 

SA  IL/3  am e>i.3 

: . nao 

INC 

57.99 

FAIL/5.4'*°L£ 

9 . *5  “.2. 

7CT 

AL  NC.  0=  SABLES  P = ? 

IN  CL  . 

Figure  17  - Sample  Output  of  the  Grader-1  Program. 


Another  feature  of  the  Grader  program  provided  t<_>  save 
computer  time  decides  which  satellites  can  be  taken  out  of 
the  dot  product  part  of  the  program.  If  one  is  analyzing 
half  of  the  northern  hemisphere,  then  there  will  be  satel- 
lites which  will  never  be  visible  to  any  of  the  position 
vectors.  A quick  survey  at  the  beginning  of  the  program  is 
thus  made  to  see  if  any  satellites  have  too  great  a negative 
'y'  component  or  negative  'z'  component  thus  making  them 
ineffective  to  the  area  in  question*.  This  is  discussed  in 
detail  in  Appendix  A. 

* Suggested  by  Robert  White  - C.S. Draper  Laboratory 


Appendix  A 


The  'Grader'  Program  in  Detail 

The  'Grader'  program  was  written  in  the  MAC  language 
which  lends  itself  well  to  vector  operations.  However,  con- 
version to  a more  common  language  should  be  straightforward. 

The  input  is  shown  in  Figure  A-l.  It  consists  of: 

'longdelt'  values  (seldom  change) 

( initial 

inclination  values  ' increment 

I final 

N — number  of  satellites 

ALPHA  — threshold  angle  of  visibility  (deg) 

F — phase  angle  number  ( 0 , 1 , 2 , 3 , . . . , N - 1) 

EA  — elapsed  angle  (deg) 

F — 

EA  — 

LONGDELT  values  are  the  incremental  arcs  in  longitude 
through  which  the  vehicle  position  vector  is  stepped.  For 
example  at  88°  latitude,  the  increment  in  longitude  for  a 
net  2°  motion  is  1.0479  radians  (or  60.04  degrees).  At  86° 
latitude,  .50588  radians  (or  28.98  degrees)  is  the  proper 
longitudinal  step  size  and  so  on.  A latitude  of  90°  does 
not  have  a corresponding  step  size  since  the  calculations 
need  only  be  done  once.  This  is  taken  care  of  by  the  pro- 
gram. How  the  numbers  were  calculated  will  be  explained 


a) 


1 

f b) 

one  J 
card  \ 

1 

|C) 

1 

<d) 

one 

(e) 

card 

if) 

one 

(g) 

card 

ih) 

later. 
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* COLLECT  MLS1 580 .GRADER 

* HE APING  NETWORK  24/24/14  AND  24/24/5. 

* GO 

R ************************************************* 

R * LAST  REVISION  TO  PROGRAM:  MAY  13.  1977 

R ************************************************* 

R LONGDELT  VALUES  (45) 

1 .0479776729  1 

5.058 859  09053E— 01 

3. 355 538 169 57E- 01 

2.51  4 94 05586 7 E—0 1 

2.01367776409E-01  Note  1 

1 .6809348852  OE  — 0 1 
1.4441S920C71E-01 
1 .26724659301E— 01 

1.1301 96831 90E-01 

1 • 0 2 1 0 32  636 1 6E  — 0 1 

9.3214249041 2E -02 

8.  584 5 73  7371 2E— 02 
7.964710 13844E-02 
7. 43680877799E— 02 
6.982531 5093 8E- 02 
6. 58  81 353  7 02  IE -02 

6.2431  1991331 E—  0 2 

5 .93932596  73  IE -0  2 

5. 67031931 536E-02 

5. 43 0959 46 74 2E- 02 
5.217  09 1 77248E— 02 
5. 0 253235954 2E- 02 
4.85285895634E— 02 
4.  69  7374  5755  OE  — 02 

4.55692571 6 75 E—0  2 

4. 429873821 5 OE -02 

4.3148302701 5E- 02 

4. 2 1 06 12 2504 3E -02 
4. 1 1620781 126E-02 
4.030747964 1 7E-02 
3.  9534 84 2462 4E-C2 
3. 8 83 7 70 5565 7E- 02 

3.821  04  836699E  — 02 

3.764834620S3E  — 02 

3.714711 7895 7E -02 

3.  67031968401E  — 02 
3.63134869029E— 02 
3.597534 19082E-02 
3.568651 96697E-02 
3. 54451 443050E-02 
3.524967561 1 7E-02 
3. 569888454 93E-6 2 

3. 4991 83  4 081 6E-C2 

3. 492 786 48 1 8 6E- 02 
3. 4 906585 03 9 8E -02 

R INC  IN, INCDELT.  INCF1N.  N,  ALPHA 

51  . 1 , 55,  24.  45 

R F.EA 

14,78 
14.80 

14,82 

14,84 

14. e6 

14.88 

lAi9  0 

14,92 
0 ,0 

40.  5*  70.  24.  45 

5.0 

Figure  A-l  - Listing  of  input  for  the  Grader  program. 
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INCIH,  INCDELT,  INCFIN  represent  the  initial  inclination, 
inclination  increment,  and  final  inclination  respectively. 
Suggested  values  are:  40,  5,  70  (coarse  plotting) 

and  50,  1,  60  (fine  plotting) . 

F,  EA,  representing  F value  or  phase  angle  value  and  the 
elapsed  angle,  are  in  the  inner  most  loop  with  a read  state- 
ment at  the  beginning  of  the  loop.  One  can  have  any  combina- 
tion whatsoever.  The  elapsed  angle  (EA)  is  the  angular  motion 
of  all  the  satellites  in  their  respective  orbital  planes. 

Alpha 

Alpha  is  the  threshold  angle.  This  earth  central  angle 
between  user  and  satellite  defines  visibility  for  the  program. 
If  a vehicle  position  vector  is  found  at  any  angle  greater 
than  this  amount  it  is  rendered  nonvisible.  For  24/24/14, 

45°  was  found  to  be  a good  number. 

It  is  important  to  choose  an  alpha  value  that  is  less 
than  the  expected  performance  of  the  constellation  so  that 
the  percentage  of  failed  samples  will  not  be  zero.  It  is 
equally  important  to  use  this  threshold  value  for  all  the 
runs  for  a given  N value  so  that  a meaningful  comparison  can 
be  made. 

More  on  Longdelt 

The  value  of  this  2°  net  longitudinal  movement,  longdelt, 
was  put  in  to  compensate  for  the  fact  that  a degree  of  longi- 
tude at  the  equator  and  a degree  of  longitude  near  the  north 
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pole  represent  radically  different  angular  values  on  the  sur- 
face of  the  globe.  If  one  were  to  disregard  this,  one  would 
find  that  a disproportionate  number  of  samples  would  be  taken 
near  the  north  polar  area. 

Let  us  represent  the  vehicle  position  vector  in  rectan- 
gular coordinates.  To  do  this  a transformation  vector  is 
needed  to  transform  from  i,  L to  x,  y,  z.  This  then  is 


cos  2,q  • cos  Lq 
cos  • sin  Lq 
sin  £q 


where  LQ  = 0. 


Then  another  position  vector,  two  degrees  away  from 

a 

pn  can  be  defined  as 


cos  Jl^  cos  L1 
cos  i sin  L^ 
sin 


where 


Z1  = *0 
L!  * ? 


Since  they  are  both  unit  vectors 


cos  9 ■ Pq  • p1 


' ' mm  **  **  *»_,•» 
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or  cos  2* 


cos  ZQ  cos  Lq  cos  Z^  cos 


+ cos  ZQ  sin  Lq  cos  Z^  sin 


+ sin  Zq  sin  Z^ 


Substituting  ZQ  for  Z^  and  0 for  Lg  yields: 


cos  2°  =»  cos  Zg  cos  Zg  cos  + sin  Zg  sin  Zg 


so  that  cos  L, 


2 

cos  2°  - sin  Z, 


cos  Z, 


L1  = cos 


For  80°  latitude  (ZQ) , 


= cos 


_1  cos  2°  - sin  Zg 


cos  Z, 


2 

cos  2°  - sin  80° 
cos2  80° 


= 11.54°  (.2014  radians.  Note  1 in 

Figure  A-l) . 


Thus  at  80°  latitude,  the  vehicle  position  vector  is  placed 
at  11.54  degree  increments  starting  with  0°,  11.54°,  23.08,. 
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At  each  latitude  then,  a corresponding  incremental  step  in 
longitude  is  made.  This  gives  a close  enough  approximation 
to  the  desired  homogeneous  sampling  pattern. 

OUTPUT 

The  output  is  straightforward,  a sample  of  which  is 
shown  again  below  for  convenience.  The  N/P/F  along  with 
elapsed  angle  describe  the  constellation  uniquely  once  an 
inclination  is  chosen.  The  'threshold  angle'  is  simply  alpha 
as  defined  above.  The  'total  number  of  samples  per  inclina- 
tion point'  is  the  number  used  to  calculate  the  'FAIL/SAMPLE' 
number.  It  is  not  computed  in  each  run  (since  it  remains  the 
same  for  this  grid)  but  was  found  once  by  means  of  a counter 
in  the  program.  The  total  number  of  failures  is  divided  by 
the  total  number  of  samples  to  get  F/S. 


N £ t * o 


2 1/2-/1  a 


C - 1 * *.  u a r;  t ?.  ?.  a * v t a « r.-  » r •*  :;  j x 

P * LAiT  -£ViiiC.\  TO  ■'.•v  12,  1?7  7 

m i t « i-  t r*.  a * :t  v t *<  « x j r t • ; « u 7»  4*r.  savnn  ^ **r  « » -.  t 
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Disregarding  Satellites  Unreachable  to  the  Sampling  Area 


The  procedure  of  disregarding  satellites  which  cannot 
be  'seen'  by  the  sampling  position  vectors  is  simple  and,  as 
mentioned  in  the  main  test,  can  save  computer  time.  What  is 
done  here  is  to  look  at  the  components  of  the  satellite  posi- 
tion vector  and  see  if  either  the  y component  or  the  z com- 
ponent is  sufficiently  negative  to  render  it  non-visible. 
Assuming  that  the  sampling  area  is  the  one  used  in  this  study, 
then 

a)  if  RRj  < sin (a),  | the  satellite  is  not  visible 

i to  the  sampling  area. 

b)  or  if  RR^  < - sin  (a)  ) 


where  RR^ 

RRk 

a 


satellite  y component 
satellite  z component 
visibility  threshold  angle  setting. 


See  notes  1 and  2 of  the  program/  Figure  A-2. 


Transformation  of  the  Satellite  Position  Vector  from 
Keplerian  Elements  to  Rectangular  Coordinates 


In  order  to  see  whether  a satellite  is  within  a specific 
angle  to  a sample  point,  the  Keplerian  elements  are  trans- 
formed to  rectangular  coordinates  on  a unit  sphere.  The  sam- 
ple point,  which  is  in  latitude,  longitude  coordinates  is 
then  also  transformed  to  rectangular  coordinates.  The  scalar 
product  of  the  two  can  then  be  taken  and  the  angle  checked. 


42 


♦OOOOIO  PMOOMAC  MLS  1 SHO. GWAOKR 
MOOOO^O  IN UK  X A.I.J.K, 
M00002S  D I MK NS  ION  (CO. 


► 


* 4 

I-  ******  * 

* * 

* 

i 

» 

* 

x — 

» 

* 

— o 

i 

* 

♦ 

u z 

X 

* 

# 

x a 

X 

# 

* 

a x 

X 2 N 

• 

* 

<\J  * 

t/5 

< c 

o 

* 

• ♦ 

z 

11.  > 

z 

* 

2 * 

< — 

• 

* 

* 

uj  a x 

X 

* 

w * 

a x 

Z X 

a 

* 

x c • 

XX  z 

2 

♦ 

— c * 

> z a 

• X - 

< 

* 

iT  * * 

C — «'il 

IT  C X 

X 

♦ 

C4  • * 

Ox 

x z a »a 

* 

• -J  * 

’SI  Z • X SI 

-<  x — c 

• c 

♦ 

VO 

D</1* 

XI  CNXtl 

CC  X V3X 

- o 

* 

0) 

1 1 

LT)  w* 

X -Wl/I  J 

2 X XU 

CM  X 

* 

c • ♦ 

— x a u o 

ox  x ax 

a 

♦ 

4* 

<J  — ♦ 

X U<J3 

x < — > 

• z 

* 

o 

x-» 

* 

X JJZ'Jin 

— — — 

* 

z 

L.  Z'/>  0 

u r - x x z 

a x 

* 

• * 

x — t/i  o 

Cx  . -3 

* X 0 

* 

-1-1* 

< a > o 

s>oz  • O - 

OX  • X 

* 

CW* 

SI  • 3 0 — 

O X<  -zx 

a - — x 

* 

• w * 

xx  C 

x z ax  j o- 

— z u < 

* 

D •* 

z z<  = o 

2 X — x vu  i/> 

■X  X X 

* 

r.  -*♦ 

Z X z c 

IX  X 3 si  C 

x > 0 2 

* 

CO* 

> 2 — X C 

a si  — x a 

w >■» 

* 

wC  ^ 

O LZ<  O 

Z x • f)  — 

X 3 • X 

-«■ 

~ • * 

iZx®? 

0 £ X U 'Ai 

- z «a  - 

•* 

c: 

* 

C x x 3 

C xX  ZX 

z < -o 

•it 

-*  w ♦ 

X Z-x  • • 

2 SI  — — 

_ X C X 

* 

< 

• ^ - 
••  -J  ♦ 

x -ckic: 

«r  r x • x x 

> • X1  Xxx 

■* 

-j 

uJ  * 

o zs'iu 

XXCDXCXx 

— x a x 

♦ 

c 

w 

a ncui 

30XX0CXU 

x C — • — X — 

«♦ 

z 

d-3* 

a ••  a x — — 

X Z < < Z -X 

z x x — z z 

* 

M 

vJ24 

ii  X 3 X X O'. 

2«  C>£  X < 

X c X O X X X 

* 

1 

> A 

w* 

CC4 

a r - z - x X 

— z x vs  x in 

— X X z > 

w> 

^ J * 

c --  a: 

X X - < > Z X 

X X C X 2 

II 

X X x ^ 

< > C — Xx 

c x < —2  a U 

♦ 

• »4 

< O X U 2 X 'X 

S,22CXQ<0 

— «fl  X X z 

* 

* 

o <zcaa 

X x XX  Z « 

20  ZE  X 

• * 

II 

-1  ‘ 

.0  rp  * 

z»x-iaa 

2 Z X x (/! 

0 2 Z Z X Z 

^5  ♦ 

X * 

x r xxx 

Oxu'xxu<x 

»cuMZCin  j* 

• IT  * 

> x (/IX  — — 

x u — x a w z 

X XX  — — z — ■*  — 

o 

• 

: <\i  O N * 

u • * 
cz  ♦ 
uc« 

■— — * 
* 

Z2» 
c X * 

— — * 


U ^IZOJ 
Li  U. X X Z 2 

uj  u o x x j 

:z-l- 


X Ll  3 < XX  ’ 


X ~ 

> IX 
X‘  X 
«r  < 2 
2 Z X' 

• < — 


• X 3 3 

■ <\J(V  X 2 s 

: x — — 

: c < z x 

. z ui  — < 

<f  2 2 
1 X 

J X II  II 


: < X x 
: xx  2 Xi 
i -xx; 


2 >CX 

zz  < 

J * 

OS  ■»  2 

c 

< U 

-c  z 

x x * 

♦ L^X 

V- 

2 Z Z 3 

_J  « 

w Ij  > 

X C< 

j 

■—  w 'w 

U»-  «J 

si  X * 

X - .X 

u 

« u 

Z < X 

O < •* 

— o ^ 

XX  x j 

-z  z 

a x * 

r>  J) 

c 

<(/'<•? 

*»• 

Z 2 S 

z 

-C  Z Z 

XX 

xu.* 

* X X 

u 

— X XI  X 

! * x t_ 
. < _x  3 
: 2 _,-(/? 


- X Q 
O < 3 
Z 3 X 
C X V- 

xzi 


X2“ <<XXX3*C  II  _ 

2 x — z < * it  c ii 

L.'  'JU'JX-Z  « i*  w y 3 

zxzxz  -c:*  <r  x r< 

— X — X — ||  .ZZ*x  c<x  x 

* X C X X £ 


* 2 

Z X 

II 

II  II  II  II  II  II  II  II  X II  II  II  * 

* X 1 

1 -< 

1 

X z * 

* a 

2 2 

< 

< < Z X — * 

* < 

c o 

r 

O X — o X X * 

* 2 

z z 

/*» 

x x x x x w jj  * r 

* U 

c c 

-i 

2 2 ZIZZZi  * - 

* 

X X 

< 

0 0 2 a-^----NZ*  c 

Q 

* 

* X 

z 

•*  ■»  * 

******* 

♦ •* 

**>*************3 

u- 

COOOoOCC  — (VOsT-TCXX 

/*  x A A <N  A ..A  A A A. 


WWW  w — v ^ X X.  — '—■WWW  WWW  WWW  w w w w W W V—W  V—  WWW  W W W w 1—J  W 

CCOOCCCOCC'COCCOOOOOOC'OOCOOOOCCOOOOOOCO 

iz:c:i3;is:c::iiJ22i2ii:c£i!i2J:rjnr'2isi."Sifl 


- 43  - 


Figure  A-2  - Listing  of  the  Grader  program  for  cases  where 
N=P.  "R”  on  the  far  left  represents  a comment 
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Figure  A-2  (continued)  The  Grader  program. 
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The  transformation  of  the  satellite  elements  then  can 
be  done  by  using  a simplified  version  of  the  R matrix  shown 


below* 


a 

-x 

*C 

qy 

= R 

qn 

_V 

-qC. 

[j  * 

M 


:■ 


I 


! 


The  figure  shown  shows  the  vector  components.  The  R matrix 
is : 


= cos  n cos  to  - sin  ft  sin  co  cos  i 
5. 2 = “ cos  ^ sin  w “ sin  ^ cos  00  cos  i 

2,^  = sin  ft  sin  i 

m^  = sin  ft  cos  go  + cos  ft  sin  go  cos  i 

m2  = - sin  ft  sin  go  + cos  ft  cos  co  cos  i 
m^  = - cos  ft  sin  i 
n^  = sin  co  sin  i 

n2  = cos  to  sin  i 

n^  = cos  i 


* Battin,  Astronautical  Guidance,  p.  18 
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— 


— 


However,  since  one  is  considering  simply  a position  on 
the  unit  sphere,  the  vector  components  and  can  be  set 
to  zero  and  q^  set  of  1.  Thus 


Notes  3,  4 and  5 on  the  program  source  show  how  RR^ , RR^ , 
and  RR^  are  derived  corresponding  to  j,^,  m^,  and  n^  respec- 
tively. 


Transformation  of  the  User  Position  Vector  to 
Rectangular  Coordinates 

Since  latitude  and  longitude  angles  lend  themselves 
readily  to  a grid  sampling  pattern  on  a unit  sphere,  they 
are  employed  here.  After  the  corresponding  increment  is 
made,  they  are  transformed  to  rectangular  coordinates. 

The  transformation  is  then. 


'V 

m 

X 

cos  (lat)  cos  (long) 

Pi 

= 

y 

= 

cos  (lat)  sin  (long) 

P2 

z 

sin  (lat) 

- - 

• m 

This  corresponds  to  notes  6,  7,  and  8 of  the  program  listing 


48 


Other  features  which  can  save  computer  time  have  been 


incorporated  and  will  be  discussed  in  Appendix  A.  By  in 


corporating  all  these,  the  author  was  able  to  cut  the  cost 


per  run  by  more  than  60%. 

Figure  A-3  shows  a listing  of  the  Grader-1  program. 
Only  the  part  which  is  different  is  shown  for  brevity.  The 
rest  is  exactly  the  same  for  both  programs. 
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Figure  A-3  - Listing  of  the  Grader-1  program  for  constellations  w’lei  e 
N>P.  Only  the  first  part  is  listed.  The  remainder  is  the 
same  as  is  found  in  Grader. 


Pigure  A-3  (continued)  The  Grader-1  program 
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Appendix  B 

Network  Orbit  Plotting  Program  and  Derivation 

After  one  has  an  idea  of  a suitable  inclination  for  a 
given  configuration  of  N satellites,  the  question  of  plotting 
the  network  and  finding  the  point  of  worst  coverage  arises. 
After  looking  at  different  ways  of  representing  the  different 
ways  of  representing  the  different  orbits  and  the  satellites 
within  those  orbits,  this  method  was  found  to  be  the  best  to 
obtain  the  quickest  spherical  plot.  From  the  plot  one  can 
then  find  the  trouble  areas  and  the  resulting  maximum  radius 
(Rmax)  for  tr^Ple  coverage. 

The  program  deals  with  both  types  of  constellations 
(studied  in  this  paper) , the  type  where  N = P and  the  type 
where  N > P,  where  P is  the  number  of  orbits.  The  N > P 
case  will  be  used  as  an  example  here.  Basically  there  are 
four  steps  to  the  program: 

1)  Input  describing  the  constellation  and  run, 

2)  the  computation  of  LAN  and  node  for  each  satellite, 

3)  the  transformation  to  rectangular  coordinates, 

4)  and  the  transformation  to  latitude  and  longitude. 

Figure  B-l  shows  the  source  programs.  It  was  developed 
on  a Hewlett-Packard  9820  desk  calculator.  Figure  B-2  shows 
the  output  of  a run  used  in  the  analysis  of  the  24/3/1  con- 
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FXD  2;  ENT  " MONTH . DATE ", R8 1 

PRT  " MONTH. DATE" ,R81;  PRT  " 1977" 

FXD  0;  ENT  "SATS" ,A 
ENT  "PLANES" ,R3 

ENT  "PHASE  NO. (F) " , R6 ; R6*360/a-  R0 

PRT  " N/P/F " , A,  R3  , R6 

FXD  1;  ENT  " INCLINATION" , B ; PRT  "INCL.",B 

ENT  "NO. OF  PTS" , R20 ; PRT  "NO. OF  PTS" ,R20 

ENT  "PLOT  ANGLE  INCR.",R21;  PRT  "PLOT  ANGLE  INCR. " 

3 6 0/A— C ? IF  A >R3 ; 3 6 0/R3  — C ; a/R3  — R7 

COS  B— Rl ; SIN  B— R2  ; 0— R4— RIO— RllO 

"10";  SPC  1;  PRT  "SAT. NO. ",R10;0— R100 

IF  A>R3 ; GTO  "13" 

R10*C— Rll ; Rl 0 * R0  —Rl 2 ; GTO  "15" 

"13";  R4*C— Rll;  R4+R0*360/R7— R12 
"15";  R12-360— R13;  IF  0>R13;GTO  "17" 

R13  — R12;  GTO  "15: 

"17";  PRT  "LAN ,NODE , Rll , R12 

Prt  " LONG,LAT " ' 

"20"; COS  Rll*COS  R12-SIN  Rll*SIN  R12*R1— R31 
IF  R31=0 ; .00001— R31 

SIN  R11*CIS  R12+COS  Rll*SIN  Rl2*Rl—  R32 

SIN  R12*R2— R-33  

ATN  (R32/R31) — R41 
IF  0 >R31 ; R41+180—  R41 
IF  0>R41;  360+R41— R41 
PRT  R41 ; PRT  ASN  R33 
R12+R21— R12;  R100+1— R100 
IF  R20>R100;  GTO  "20” 

R110+1— R110;  IF  R110>R7-1;  0— R110;  R4+1— R4 
R10+1— R10;IF5R10  A-l;  GTO  "10" 

SPC  7;  END 


1 


2 


3 


4 


Figure  3-1  - Program  giving  latitude  and  longitude  in  time  of 

multisatellite  constellations. 
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Figure  B-2  Part  of  the  output  for  plot  of  24/3/1  on  a globe. 
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stellation.  Let  us  first  look  at  the  output  closely  to 
understand  the  purpose  of  the  program.  Then  the  algorithm 
which  gives  that  output  will  be  discussed. 

The  Output 

Month-Date  - 5.02  - 1977:  the  program's  first  question 
to  the  user  is  the  date.  The  number  5.02  was  typed  in  refer- 
ring to  May  2. 

N/P/F  - 24,3,1:  This  refers  to  the  constellation  24/3/1. 

INCL.  - 53.0:  The  inclination  is  requested  by  the  pro- 
gram and  then  printed. 

NO.  OF  PTS.  - 5.0:  The  number  of  points  to  be  plotted 
for  each  satellite  can  be  expressed  also  as  the  number  of 
nodal  positions  a satellite  will  inhabit  in  its  orbit. 

PLOT  ANGLE  INCR.  - 10.0:  The  plot  angle  increment 
refers  to  the  increment  in  nodal  angle  for  these  positions. 

In  this  case,  therefore,  we  are  looking  at  node  angles  of 
0,10,20,30,  and  40°. 

SAT.  NO.  - 0.0:  This  refers  to  the  satellite  number. 

Zero  is  used  because  it  makes  the  programming  simpler. 

LAN,  NODE  - 0.0  - 0.0:  These  two  numbers  represent 
the  longitude  of  ascending  node  and  the  node  angle  of  the 
n-th  satellite. 

LONG,  LAT  - 0.0,  0.0;  6.1,  8.0;...:  This  is  the  desired 
output  and  it  is  the  longitude  and  latitude  of  the  satellite 
at  the  node  angles  requested.  See  note  B-l. 
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One 

is  then  ready 

to  plot. 

The 

following  is  a 

description  of  the  variables  found 

in  the  program. 

HP 

Notation 

Regular 

Notation 

Description 

A 

N 

the  number  of  satellites 

B 

i 

inclination  (this  assumed  the 
same  for  all  the  planes) 

C 

1 

An 

increment  in  longitude  g ascend- 
ing node  from  plane  to  plane 

RO 

Aw 

increment  in  node  angle 
(Aw  = F x 360/N) 

R1 

cos  i 

R2 

sin  i 

R3 

P 

number  of  planes  in  which  the 
satellites  are  distributed 

R4 

the  number  of  individual  plane 
(0,1, . . . , P - 1) 

R6 

F 

the  phase  angle  number 

R7 

NPO 

number  of  satellites  per  orbit 

NPO  = N/0 

R8 

<J> 

phase  angle  between  a given  satel- 
lite and  the  satellite  in  the  next 
adjacent  longitude  of  ascending 
node,  <J>  = F x 360/0 

RIO 

i 

satellite  number 

Rll 

fli 

longitude  of  ascending  node  of 
ith  satellite 

R12 

ui 

node  of  the  ith  satellite 
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HP 

Notation 

Regular 

Notation 

Description 

R1 

.3 

dummy  variable  for 

R31 

R32 

R33 

R41 

R81 

P.110 


per  satellite 

increment  through  which  all  the 
satellites  are  to  be  moved  in 
their  respective  orbits 

x component  of  the  ith  satellite 

y component  of  the  ith  satellite 

z component  of  the  ith  satellite 

longitude  of  the  ith  satellite 

month. date  in  one  variable 

jth  satellite  of  a given  orbital 
plane  when  N > P (0,1,2,..,  NPO  - 1) 


pmmmm  wmmm 


Appendix  C 


Program  for  Plotting  R vs.  Inclination 
For  a Given  Group  of  Satellites 


Once  the  satellites  surrounding  the  worst  coverage 
points  have  been  identified  as  in  Section  III,  the  program 
for  plotting  R versus  inclination  is  very  helpful.  This  was 
shown  in  Section  III  with  the  example  of  the  constellation 
24/24/14  where  it  was  used  to  produce  Figures  13  and  14. 
Although  it  was  not  used  in  obtaining  the  final  results 
given  in  Table  1 (p.  9),  for  better  accuracy  something  sim- 
ilar might  be  employed. 

Listed  below  are  the  steps  in  the  program  originally 
written  for  a large  desk  calculator  (Hewlett-Packard  9820) . 


w.. 

4: 

ENT  "LflNi NODE-Sfl 

v. 

T B ” 5 R22  .•  R42  5 

r'r.l 

...  i r — 

FRT  ” SRT  B " ? R 2 2 ? 

"h 

R42h 

1 : 

| 

3! 

FXD.  0 

SENT  "SPIT  N 

| ENT  "LAN* NODE  Sfi 

OS?"  j 

R8  t R?  * R 1 0 > 

T C " ■ R 2 3 j R 4 3 5 

FRT  " 

S ft T NOS 

FRT  "SRT  C " « R 2 3 ? 

a a a a 

? R 8 « R 9 j P.  1 0 1-  ! 

R43H 

hi  • 

6: 

FtfB  2 

> ENT  "MONTH 

ENT  "INITIAL  INC 

. BATE 

? " > R 6 5 P R T " 

?m»b;ent  “inc.de 

DATE " j R6 5 

LT7"jC?ENT  "FI  nr 

FRT  " 

L INC?" ? R 7 h 

1977" 

h 

^ a 
( • 

3 • 

FRT  " IN.  - BELT-  S.F 

ENT  " 

LAN»N0DE?  S 

INRL  I NC  " ? F'RT  "I 

FIT  R" 

> R21 » R41 ? 

NCLI NATION" « B - C - 

FRT  " 

LAN- NODE-Sfl 

R7h 

T R"  , 

R 2 1 j R 4 1 h 
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E N T "EL  R P S E E N 0 L 
E?” < R 1 3 5 FRT  "ELR 
FEED  NODE"  ? P 1 3H 
3 ; 

O+RSSPC  1 5 F RT  "L 

0 N G S « LRTS"  5 B + Rlh 
10! 

"20"  5 P i'21+fi)  + R1  1 ‘ 

5 P ( 4 1 + fi  l + R 1 3 •*  R l 2 
5 COS  Rll+COS  R 1 2 
-SIN  R 11  SIN  R12* 
COS  Rl+R (31+fl) F 
1 1 : 

IF  R (31+fl) =05 1E- 
9 + R (31+fl) V 
12:  i 

SIN  R 1 1 * C 0 S R12+  i 
COS  R11*SIN  R12* 
COS  R 1 - R ( 7 1 + fl ) h 
13* 

SIN  R 1 2 * S I N R 1 +R 

1 9 1 + fl ) P j 

14: 

RTN  (R (71+fl) ,'R (3 
1 + R ) J 4 R ( 3 1 + fl ) h 

1 *1  ■ 

IF  0>R  (31+fl) 5 R (5  ! 

1+fl) +180*R(51+fl) 
h I 

16: 

IF  0>R  (51+fl)  5 360 
+ R (31+fl) +R  (51+fl) 

■ ' i 

PRT  R (51+R) 5 RSH  • 1 
R (91+fl) +R (61+fl) 5 
FRT  R (61+fl) H 
; is: 

R+ 1 +fl  5 I F R£25  I 

GTO  " 20  " R 


1?: 


2 0 * . 

HCS  (P32P33+R72R 
73+R92R93) + R2h 

21  : 

RCS  ( R3 1 R33  + R7 1 R 
73+R91R93) + R14P 
22: 

RCS  (R31R32+R71R 
72+R91R92) + R15F 
23: 

RCS  ((COS  R2- 
COS  R14+COS  R 1 5 ) 
/(SIN  R14+SIN  R 1 
5) ) +R3H 
24: 

RCS  ((COS  Rl^- 
COS  R 2 * C 0 S R 1 5 ) / 
(SIN  R 2 * S I N R 1 5 ) 

J +R4(- 
23: 

RCS  ((COS  R 15- 
COS  R2+C0S  R 1 4 ) / 
(SIN  R 2 * S I N R 1 4 ) 

) + R3h 
2S  • 

RTN  (TRN  ( R2/2 ) * 
1/COS  ( . 5 ( R4+R5- 
R3  J J ) + R 1 O 5 h 
27 : 

FRT  "FOP  I N C . = • . 

. . . . m»ri;prt  "Rfl 

DIOS  = "»R 

1 03  5 SRC  21- 
28 : 

Rl+C  + Rl 5 IF  RLLR7- 
5 0 + R 5 GTO  "20“H 
2? : 

SFC  6h 


Below  is  an  explanation  of  the  variables  found  i; 
the  program.  The  units  are  in  degrees  unless  otherwise 
specified. 
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Program  Corresponding 
Variable  Variable 


Comments 


A 

counter  (no  units) 

B 

10 

initial  inclination 

C 

Ai 

increment  in  inclination 

R1 

i 

inclination 

R2 

a 

side  a of  the  spherical 
triangle 

R3 

A 

angle  A opposite  to  side  a on 
spherical  triangle 

R4 

B 

angle  B opposite  to  side  b on 
spherical  triangle 

R5 

C 

angle  C opposite  to  side  c on 
spherical  triangle 

R6 

month. date  of  the  run 

R7 

final  inclination 

Rll 

°i 

longitude  of  ascending  node 
for  the  i-th  satellite 

R12 

OJ  . 

l 

node  angle  of  the  i-th  satellite 

R13 

EA 

elapsed  node  under  investigation 

R14 

b 

side  b of  the  spherical  triangle 

R15 

c 

side  c of  the  spherical  triangle 

R21 

fiA 

longitude  of  ascending  node  of 
satellite  A 

R22 

°B 

longitude  of  ascending  node  of 
satellite  B 

R23 

nc 

longitude  of  ascending  node  of 
satellite  C 

R31 

XA 

component  in  the  x direction  of 
satellite  A (length  units) 

Program 

Variable 


Corresponding 

Variable 


Comments 


«i 

i 


R32 

XB 

component  in  the  x direction 
satellite  B (length  units) 

of 

R33 

xc 

component  in  the  x direction 
of  satellite  C (length  units) 

of 

R41 

“A 

node  angle  of  satellite  A at 

0 E.A 

R42 

“b 

node  angle  of  satellite  B at 

0 E.A 

R43 

wc 

node  angle  of  satellite  C at 

0 E.A 

R51 

*A 

longitude  of  satellite  A 

R52 

*B 

longitude  of  satellite  B 

R53 

*c 

longitude  of  satellite  C 

R61 

la 

latitude  of  satellite  A 

R62 

latitude  of  satellite  B 

R63 

Lc 

latitude  of  satellite  C 

R71 

yA 

component  in  the  y direction 
satellite  A (length  units) 

of 

R72 

yB 

component  in  the  y direction 
satellite  B (length  units) 

of 

R73 

yc 

component  in  the  y direction 
satellite  C (length  units) 

of 

R91 

ZA 

component  in  the  z direction 
satellite  A (length  units) 

of 

R92 

ZB 

component  in  the  z direction 
satellite  B (length  units) 

of 

R93 

zc 

component  in  the  z direction 
satellite  C (length  units) 

of 

RIO  5 

R 

the  final  resulting  radius  of  the 
circumcircle  containing  satellites 
A,  B,  and  C on  its  circumference 
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The  Input 


The  input  to  the  program  is  shown  as  typed  by  the  com- 
puter in  Figure  C-2. 


1 


I 


OUTPUT 


STP 

SAT  NOS 

u 

4 

19 

PATE 

4.05 
1 977 

LAN  j NODE-SAT  A 

30.00 

60. 00 

SAT  B 

6 0 . 0 0 
1 20. 00 

SAT  C 

285.00 

30.00 

IN.  » PELT j &F INAL 
INCLINATION 

50.00 
1 . 00 

60. 00 

ELAFSED  NODE 

15.00 


LONGS  > L AT 


FUR  I N C . = 
RAH  I US  = 


97 

47 


80 


50 . 00 


FUR  I N C . = . 
RADIUS  = . 


96. 94 
48.65 

fcl'  0 f • 8 


’ 51  ! 00 


51.24 


FOR  INC. = 
RADIUS  = 


96. 48 
49.57 

o o o 
£_  u 

•o  o o *- 
•»'  •-.«  • w O 

316.62 
3 3 . 8 6 

52. 08 

■ ■ • • • 

50. 53 


96. 00 
50.48 


208.96 
34.3  8 

316.04 

34.38 

FUR  INC.=, 

53.00 

F.RDIUS  = , 

4 9.83 

Figure  C-2  - Sample  Input  and  Output  of  the  R vs.  Inclina- 
tion Program.  It  gives  the  radius  (R)  of  the 
circumcircle  containing  satellites  2,  4,  and 
19  on  its  circumference  for  varying  inclina- 
tions. .It  also  gives  the  corre'fepong  longitudes 
and  latitudes.  See  graph  on  Figure  13. 
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RRDIU 


RRDIU 


RRDIUS 


RFID  I US 


RRDIUS 


RRDIUS 


Figure  C-2  (continued) 


Input  #1  (2,4,19)  - the  numbers  of  the  satellites  comprising 

this  spherical  triangle  as  taken  from 
Appendix  B.  (only  for  reference  purposes) . 

Input  #2  (4.05)  - this  is  the  month  and  the  date  of  the  run. 

Input  #3  (30,60)  - the  longitude  of  ascending  node  (Q)  and 

node  angle  (u)  of  satellite  A as  found 
in  Appendix  B. 


Input  #4  (60,120)  - the  same  as  input  #3  but  for  satellite  B. 

Input  #5  (285,30)  - the  same  as  input  #3  but  for  satellite  C. 

Input  #6  (50,1,60)  - these  are  the  initial  inclinations  of  the 

satellite  orbits,  the  increment,  and  the 
final  inclination. 


Input  #7  (15)  - the  elapsed  node  of  all  the  satellites  in 

their  orbits. 


The  output 

A sample  output  is  listed  above  taken  from  the  24/24/14 
constellation.  The  satellites  involved  are  2,  4,  19,  21  (9) 
(14)  as  listed  in  Figure  14.  The  three  used  for  the  program 
are  2,  4,  19.  Only  three  of  the  four  satellites  are  required 
to  define  the  circle  whose  radius  is  being  sought.  One  must 
be  sure  that  the  fact  that  there  were  four  satellites  equi- 
distant from  a given  point  is  not  coincidental  only  to  the 
particular  inclination  at  which  it  was  analyzed.  In  other 
words  it  might  be  that  satellites  2,  4,  19  and  21  do  not  all 
remain  on  the  circumcircle  as  the  inclination  is  varied.  This 
was  checked  in  the  above  example. 
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; 


•--vm 


Essentially  lines  10  through  13  compute  the  rectangular 
coordinates  of  the  satellites  from  their  Keplerian  elements 
as  was  done  in  the  Grader  program. 

Lines  14  through  18  compute  the  longitude  and  latitudes 
of  each  satellite  at  each  inclination  for  reference. 


.ft 


Lines  20  through  22  compute  the  angular  length  of  the 
sides  which  make  up  the  spherical  triangle — sides  a,  b,  and 
c.  This  is  done  by  taking  the  dot  product  of  the  correspond- 
ing satellite  position  vectors  and  then  taking  the  arccosine 
to  find  the  earth  central  angle  between  them.  This  angle  is 
then  the  arc  comprising  one  side  of  the  spherical  angle. 

Lines  23  to  25  employ  a trigonometric  formula  to  give 
the  angles  of  the  spherical  triangle  once  the  three  sides 
are  done.  The  formula*  states  that: 


cos  A = 


cos  a - cos  b cos  c 
sin  b sin  c 


This  in  turn  is  also  used  to  produce  angles  B and  C 


Line  26  gives  the  radius  of  the  circumcircle  (R)  which 
contains  the  three  satellites  in  its  circumference.  This 

4 

handy  formula  is : 


tan  R = 


tan  -j-  • sec  -y-(B  + C - A). 


* Taken  from  standard  Math  table. 
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